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Abstract 

We consider the initial boundary value problem for the focusing nonlinear Schrodinger equa- 
tion in the quarter plane x > 0, t > in the case of decaying initial data (for t = 0, as x — > +oo) 
and the Robin boundary condition at x = 0. We revisit the approach based on the simul- 
taneous spectral analysis of the Lax pair equations [TD], [T3] and show that the method can 
be implemented without any a priori assumptions on the long-time behavior of the boundary 
values. 

1 Introduction 

The inverse scattering transform (1ST) method for studying initial value (IV) problems for certain 
nonlinear evolution equations — integrable nonlinear equations possessing a Lax pair representation 
- is known as a powerful tool for obtaining rigorous results concerning the most subtle issues 
about the behavior of solutions of these problems, including detailed long-time behavior. The most 
efficient implementation of the 1ST method turns out to be based on the Riemann-Hilbert (RH) 
problem method, which is essentially the reformulation of the scattering problem for one of the Lax 
pair equation — the x-equation — in terms of an analytic (matrix) factorization problem of the 
Riemann-Hilbert type. 

The studies on the adaptation of the 1ST method to initial boundary value (IBV) problems 
yield particular classes of boundary conditions, under which the IBV problem remains completely 
integrable, i.e. solving it reduces to solving a series of linear problems. At the beginning of 1990s 
it has been realized |17} [TBI EJ [TJ [19] that this is the case when the boundary conditions allow an 
appropriate continuation, based on the Backlund transformation, of the given initial data to the 
whole axis, which reduces the study of the IBV problem to the study of the associated IV problem. 

In a recent paper [5], this Backlund transformation approach is thoroughly revised and adjusted 
to the modern, Riemann-Hilbert framework. This allowed the authors of [5] to apply the nonlinear 
steepest descent method [6] and analyze in great details several important asymptotic and stability 
questions |16j in the theory of the focusing nonlinear Schrodinger equation (NLS). 

A general approach to IBV problems for integrable nonlinear equations was originated in 1993 
in [9], and it has been actively developed since then [HI Q21 HJ [15] . In this approach, one treats the 
IBV problems with general boundary conditions by using the spectral analysis of the both linear 
equations constituting the Lax pair. The relevant Riemann-Hilbert problem is posed now on the 
union (cross) of real and imaginary axis, rather than on the real axis only as in the case of the 
IV problems. A key factor that affects the efficiency of the method is that the construction of 
the underlying Riemann-Hilbert problem requires the knowledge of the spectral functions that are 
determined, in general, by an excessive number of boundary values. These boundary values cannot 
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be prescribed arbitrarily for a well-posed IBV problem and thus the problem of compatibility of 
the boundary and initial values arises. It turns out that this compatibility can be expressed, in 
a rather simple, explicit way, in terms of the associated spectral functions [15j . which in turn 
allows obtaining, particularly, the detailed asymptotic pictures. But translating this description of 
compatibility into the physical space (of boundary and initial values) requires solving, in general, 
nonlinear problems [3], which makes the whole problem non-integrable. 

However, for particular boundary conditions, an additional symmetry in the spectral space 
allows bypassing the nonlinear step of resolving the compatibility issue and thus making the IBV 
problem as integrable as the associated IV problem [10] (see also [H [151 03] ) • 

In [TO] (see also [15] ) , the boundary conditions (called linearizable) leading to integrable IBV for 
the NLS equation have been specified by applying a symmetry analysis to the associated Lax pair 
equations. Particularly, the Dirichlet, Neumann, and Robin boundary conditions were selected in 
this way, and the associated Riemann-Hilbert problem formalism was presented, under assumption 
that the boundary values being considered as a "potential" in the time-type equation from the Lax 
pair decay as t — > oo and, moreover, generate the associated spectral functions with appropriate 
analytic properties. Furthermore, in the paper |13] it was shown that for the integrable boundary 
conditions and, specifically, for the Robin boundary condition, the original master RH problem 
posed on the cross can be unfolded to the problem posed on the real line only. This, in principle, 
allows one to reproduce the earlier results of [19] within the general scheme of [9j. 

The important issue which has not been fully addressed in [10], [15] . and [13] is the following. 
As it has been already mentioned, the original master RH problem is formulated on the cross under 
a priory assumption that the boundary data rapidly decay as t — > oo. On the other hand, using 
simple explicit soliton-like solutions of the NLS equation, one can easily see that this assumption 
is not necessarily true for solutions satisfying a linearizable boundary condition. Indeed, the two- 
parameter solutions 

u(x,t) = — e ia * 

cosn(ax + (po) 

with a € M and </>o £ K, which are obviously not decaying, for any x, as t oo, satisfy the Robin 
condition u x (0, t) + q ■ u(0, i) = with q = a ■ tanh^o- Moreover, in the case q > the large-t 
behavior of the solution is oscillatory for generic initial data (see the discussion in Section 4). Hence 
the necessity to produce an independent proof that the Riemann-Hilbert problem obtained in |15j 
and [13] for linearizable boundary conditions indeed yields the solution of the IBV problem in 
question. In the present paper, using some ideas which go back to the late 1980s works on algebro- 
geometric solutions of integrable PDEs, we show how this independent proof can be achieved and 
hence complete the program originated in |10| . 

We revisit the IBV for the NLS equation with linearizable (Robin) boundary condition 

l 1 9 

iut + u xx + 2\u\ u = 0, x > 0, t > 0, 

u(x, 0) = uo(x), x > 0, (1.1) 

u x (0,t) + qu(0,t) = 0, t > 0, 

where (a) uq(x) decays to as x — > +oo and (b) q is a real constant. Solving (II. ip by the RH 
method consists of three steps: 

(i) provide a family of the Riemann-Hilbert problems parametrized by x and t such that the 
solution u(x, t) of (11.1 j) is expressed in terms of the solutions of these problems; 
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(ii) prove that u satisfies the initial condition u(x,0) = Uq(x); 
(hi) prove that u satisfies the boundary condition u x (0,t) + qu(Q,t) = 0. 

Now let us comment on this procedure. Concerning (i), the construction of the RH problem 
has to involve only the spectral functions associated with the initial condition uo(x). In the general 
framework of the simultaneous spectral analysis of the Lax pair equations, the RH problem is 
naturally formulated (see [15]) on the contour consisting of two lines: k € M and k € iR; this reflects 
the fact that the spectrum of the t-equation from the Lax pair with coefficients that are finitely 
supported or decaying at infinity consists of these lines. But for particular boundary conditions, 
the contour can be deformed (two rays of the imaginary axis can be fold down to the positive real 
axis) to the real axis only [13]. Then it is the inherited symmetry property of the jump matrix for 
the deformed problem that allows verifying directly that the boundary condition holds. As for the 
verification of the initial condition, it is based on the fact that for x = 0, the original RH problem 
can be deformed in the opposite way (the rays of the imaginary axis are fold down to the negative 
real axis) thus reducing the problem to that associated with the x-equation of the Lax pair with 
the potential uo(x) [T5] . 

In this way we re-derive the results of |19] , [5] and hence show that the approach to linearizable 
IBV problems stemming from the general methodology of simultaneous spectral analysis of the Lax 
pair equations |10j . |15j can be implemented without any a priori assumptions about the long-time 
behavior of the boundary values. 



2 The Riemann— Hilbert formalism for IBV problems 

First, let us recall the Riemann-Hilbert formalism for IBV problems on the half-line x > for the 
NLS equation [IB] . 

The focusing NLS equation 

iiij + u xx + 2|it| 2 w = (2-1) 
is the compatibility condition of two linear equations (Lax pair) [20j: 

V x + ika^ = (2.2) 

with 

" - o) < 2 - 3 > 

and 

% + 2ik 2 a 3 ^ = (2.4) 

with 

V = ( l][U \ 2 2ku + 'm x \ , . 

^-2^ + iu x -i|n| 2 ) ^ ' ' 

Assuming that u(x,t) satisfies (|2.ip for x > and < t < T with some T < oo, define 
the solutions ^j(x,t,k), j = 1,2,3 of fl2^D- rf2~5]) as follows: Vj(x,t,k) = <f>j(x, t, k)e^~ ikx - 2ik2t ^ 3 , 
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where <3?j solve the integral equations 

t,k) = I- e~ ikxa3 J e - 2ifc2 ('- r ) <T3 y(0, r, fc)$i(0, r, fc) e 2ifc2 (*- T ) <T3 e ifcxCT3 

+ [ X e- ik ( x -y^U(y,t)^ 1 (y,t,k)e lk{x - y)a3 dy, (2.6a) 

JO 

/•t 

$2(x,t,k) = I + e~ ikx ^ / e- 2ifc2( '- r)<T3 y(0,r,/t)$2(0,r,fc)e 2ifc2(t - r),T3 e ifcx<J3 

■/ o 

+ I e- ik ( x - y ^U(y,t)<P 2 (y,t,k)e ik ( x - y ^dy, (2.6b) 
J o 

/>oo 

$ 3 (s, t,k) = I- e- ik ^ x -^U(y, t)$ 3 (y, t, k^^dy. (2.6c) 

Define the scattering matrices s(k) and >S(/c), A: £ M, as the matrices relating the eigenfunctions 
&j(x, t, k) for all a; and i: 

V 3 (x,t,k) =V 2 (x,t,k)s(k), 9i(x,t,k) =9 2 (x,t,k)S(k). (2.7) 

The symmetry 

^n(x,t,k) = ^ 22 (x,t,k), y 12 (x,t,k) = -V 21 (x,t,k) (2.8) 

implies that 

From (|2.6p and (|2.7p it follows that the spectral functions a(&) and b(k) can be analytically continued 
into the upper half-plane k £ <C + as bounded functions, with a(k) — > 1 and b(k) — > as k — > oo. 
Moreover, they are determined by it(x, 0), x > only, via 

= (^ 2 ) _1 (a;,0,fc)^3(x,0,A;). 

Similarly, the spectral functions A = A(k;T) and B = B(k;T) are entire functions bounded in 
the domains I and III, where I = {k : Imk > 0,ReA; > 0} and III = {k : Irak < 0,Re&: < 0}, 
with A(k;T) — > 1 and B(k;T) — > as k — > oo. Moreover, they are determined by u(0,t) and 
u x (0, t) for < t < T only, via 

S(k;T) = (^ 2 )- 1 (0,t,A;)^ 1 (0,i,A ; ). 

The compatibility of the set of functions {u(x,0),u(0,t),u x (0,t)} as traces of a solution u(x,t) 
of the NLS equation can be expressed in terms of the associated spectral functions as follows: 

A(k;T)b(k) -a(k)B(k;t) =c(k;T)e 4ik2T , lmfc>0, (2.10) 

with some analytic c(k;T) = O(^) as k — > oo (in the general scheme |10| of analysis of IBV 
problems, fl2 . lQf) is called the global relation). 

Define 

d(k) = a{k)A(k) + b(k)B(k), k e II = {k : Imk > 0,Rek < 0}. (2.11) 
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Finally, assuming that d(k) has at most a finite number of simple zeros in II and a(k) has at most 
a finite number of simple zeros in C + , define a piecewise meromorphic function (the superscripts 
denote the column of the respective matrix) 



M(x,t,k) 



(i) 



(2) 



(2) 



5" *?, 



(1) 



2) 



Im k > 0, Im A; 2 > 
Im fe > 0, Im k 2 < 
Im < 0, Im k 2 > 
Im < 0, Im k 2 < 



(2.12) 



Then the scattering relations (|2.7|) imply that the limiting values of M on the cross Imk 2 = 
satisfy the jump relations 



M + (x,t,k) = M-(x,t,k)e 



-i(kx+2k 2 t)a 3 



J Q (k)e 



where 



Jo(fc) 



l+|r(A;)| 2 f(Jfe) 
r(fc) 1 
1 N 

i f(k) 

1 

1 + \r(k) + T(A;)| 2 r(jfe) + T(Jfe) ' 
r(jfe) + r(jfe) 1 



i(fcx+2fe 2 i)(T3 

k > 0, 
A; € il 
fe G il 
A; < 0. 



(2.13) 



(2.14) 



where r(A;) = b(k)/a(k), 



F(k) 



B(k) 



(2.15) 



a{k)d{k) 

The orientation of the contour is chosen as from -co to +oo along R and away from along iR. 

Complemented with the normalization condition M = I + 0(l/k) as k — > oo and the respective 
residue conditions at the zeros of a(k) and d(k) (see [15] for details), the jump relation (12. 13ft can 
be viewed as the Riemann-Hilbert problem: given {a(k),b(k), A(k), B(k)}, find M(x,t,k) for all 
x > and t > 0. Then the solution of the NLS equation, u(x, t), is given in terms of M(x, t, k) by 



u(x,t) = 2i lim kM^ix, t, k). 

k— >oo 



(2.16) 



Moreover, u(x, 0) generates {o(A;), 6(/c)} and {u(0, t), u x (0, t)} generates {^4(/c), 5(A;)} as the cor- 
responding spectral functions provided the latter verify the global relation (I2,10p , Therefore, the 
Riemann-Hilbert problem approach gives the solution of the overdetermined IBV problem 



iut + u xx + 2\u\ u = 0, 
u(x,Q) = uo(x), 
u(0,t) = v (t), 
u x (0,t)=v 1 (t), 



x > 0, t > 0, 
x > 0, 
< t < T, 
< f < T 



(2.17) 
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provided that the spectral functions {a(k),b(k), A(k), B(k)} constructed from {uq(x), v (t), vi(t)} 
satisfy the global relation (|2.10p . 



3 The Riemann— Hilbert formalism for Robin boundary condition 

Even in a conditional context presented in the previous Section, the RH method allows obtain- 
ing useful information about the solution, e.g., that concerning the large-time behavior, see |15j . 
However, there are cases when one can overcome the conditional nature of the solution and solve 
a well-posed initial boundary value problem. The key factor [10j(see also [15J for making this 
possible is the existence of an additional symmetry in the spectral problem for the i-equation of 
the Lax pair. 

This holds in the case of Robin boundary condition. Indeed, if u + qu x = with some q E R, 
then the matrix V := V — 2i/c 2 o"3 of the i-equation ^ = V^> satisfies the symmetry relation [15] 



V(x,t,-k) = N(k)V(x,t,k)N- 1 (k), (3.1) 



where N(k) = di&g{Ni(k), N 2 (k)} with Ni(k) = 2k + iq and N 2 (k) = -2k + iq. In turn, ([HTTP 
implies the symmetry for S: S(—k;T) = N(k)S(k;t)N~ 1 (k), which reads in terms of A and B as 

A(-k;T) = A(k;T), 

B(-k;T) = -^±^B(k;T). (3.2) 

Now observe that the global relation (|2,10p combined with the symmetry relation (|3.2|) allows 
rewriting the RH problem ()2.13|) in the form that uses only the spectral functions a(k) and b(k) 
associated with the initial values u(x,0). Indeed, since the exponential in the r.h.s. of (|2.10p is 
rapidly decaying for k E /, the global relation (|2.10p suggests to replace -f (fc; T) by \(k) for k E I. 
Then, the symmetry relation (|3.2|) suggests to replace ^(k;T) by — |f^f( — k) for k E III and 

consequently to replace ^(k) T) by — f^zj||( — k) for k E II, including the boundaries of respectively 
III and II. The resulting jump conditions have the same form as in (|2.13p . 

M+(x,t,k) = M_(x,t,fc)e~ i(fe:i ' +2fc2 * )ff3 Jo(/fc)e i(fcx+2fc2 ' ),T3 , (3.3) 
but with T(k) replaced (cf. [15]) by 



f (k) = b{ ~ k) 2k + ' q (3 4) 

a(k) (2k - iq)a(k)a(-k) - (2k + iq)b(k)b(-k) ' 

Now we note that although T(k) is defined, for general boundary values, only for k E iT, the 
function T(k), in the generic case of absence of the zeros of the denominator in (13. 4h . is analytic 
(and bounded) for all k E C + (and continuous up to the boundary). On the other hand, the 

/ i o\ (\ f(k)e~ 2ikx ~ Aik2t \ 

exponentials in ( 2ifcic+4ifc 2 t ]) anc ^ Iq I are bounded in respectively I and 

IV. Thus we can deform the RH problem with jump (|3.3p on the cross to that on the real axis by 
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defining 

'M{x,t,k), kellulll, 



M(x,t,k) = < 



M(x,t,k) I ~ J, , A ., 2 + ° | , k G /, 

K ! \r(k)e 2lkx+4lk t 1 ) (3.5) 



Ac.,**) (J - r «*™' 1 , * 6JV . 



The resulting jump conditions take the form 

M+(x,t,k) = M_(x,t,A;)e- i(fea;+2fc2 * ),73 Jo(A;)e i(fca;+2/c2 * )ff3 , kR, (3.6) 



where 



with 



Mk)=( 1+ ^f f f T ) (3.7) 



r.(fc) = r(k) + f (*) = (2fc-ig)&(fe)a(-fe) + (2fc + W-fe)aW (3 . g) 

(2A; - i9)o(fc)o(-A;) - (2k + iq)b{-k)b(k) 

If the denominator in (|3.8|) has zeros in C + , then the formulation of the Riemann-Hilbert problem, 
normalized by M — > I as k — > oo, is to be complemented by the residue conditions at these points 
(in this case, we make the genericity assumption that these zeros are simple and finite in number). 

Remark 1. As it has already been mentioned in the introduction, the important observation that 
under the symmetry relations (|3.2p the RH problem on the cross can be deformed to the RH 
problem on the real axis was first made in |13j . 

Remark 2. We emphasize that the analytical continuation T(k) of the function T(k) is not obliged 
to satisfy equation (|2.15p on the positive real axis where, in view of the global relation, it would 
have led to the erroneous conclusion that r e (k) must vanish for all positive k. 

The RH problem (|3.6j) coincides with the RH problems obtained in |19j and [5] via the Backlund 
technique mentioned in the introduction. Our derivation, as being based on the general IBV 
methodology, is different. It also should be noted that we have made certain a priory assumptions. 
Indeed, to be able to push T — > oo in (|2.10p and hence replace ^(k;T) by |(fe) for k G / we need 
to assume a fast i-decay of the boundary data. The latter is not necessarily true for the Robin 
boundary conditions (indeed, it is generally not true). In what follows we show that, independently 
of the previous considerations, the RH problem (|3.6p yields indeed the solution of the NLS equation 
on the half line with the Robin boundary condition. To this end, as it has already been indicated in 
the Introduction, we shall show that the function u(x,t) generated by this RH problem satisfies (a) 
the NLS equation, (b) the given initial conditions, and (c) (most challenging) the Robin boundary 
condition. 

The first part of the program is easy. The RH problems (|3.3p and (j3.6|) both give the solution 
of the NLS equation in the domain x > 0, t > via (|2.16p : this is a standard fact based on ideas 
of the dressing method, see, e.g., [7]. 

In order to verify the initial condition u(x,0) = uq(x), one observes that for t = 0, the expo- 



1 ^andf 1 r(k)e- 2ikx 
T(k)e 2ikx l) and U 1 



nentials in ( ^ n ^ 2 \kx i ) anc ^ I -, I are a ^ so bounded in respectively II and III. Thus 
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we can deform the RH problem with jump (|3.3p on the cross to that on the real axis by defining 



M(x,t,k) 



r M(x,t,k), kelUlV, 
1 N 

Fe 2ikx 1) ' ~ (3.9) 



M(x,t,k) [ „ 2 . kx J , ken, 



1 fe~ 2ifc:r 1 \ 
M(x,t,k) | I , fc G III, 



which results in the jump condition 

M + (x,0,fc) = M_(x,0,A:)e~ ifc:E<T3 Jo(A:)e i(fc:E<T3 , keR, (3.10) 

where 

^)=( 1 + ll* )|2 «*>Y (3.1D 



r(fc) 1 

But the resulting RH problem (with residue conditions modified appropriately, see [15]) coincides 
with that for the spectral mapping {uo(x)} >-» {a(k),b(k)}, which yields u(x, 0) = uo(a;) due to the 
uniqueness of the solution of the RH problem. 

The usefulness of the RH problem in the form (j3.6|) is that it allows to verify that u(x, t) satisfies 
the Robin boundary condition by using a symmetry of r e (k) followed from its construction (|3,8p . 
This symmetry is as follows: 

r e (-k)=r e (k)^L (3.12) 

where 



a(k) = (2k - iq)a{k)a{-k) - (2k + iq)b(k)b(-k). (3.13) 

It is convenient to normalize a(k), which is analytic in C + , in such a way that it approaches 1 
as k — >• oo and that it has neither a zero no a pole at k = Hp. Depending on the sign of q and 

M 

2 ' 



the behavior of a(k) and 6(/c) at k = Hp, different normalizing factors are needed. Indeed, if one 



introduces a e (k) and j3 by 

( M _ JaS = ^M^) ~ j^6(fe)&Ffe), if g<0,a(-f)^0or g>0,6(f) = 
\^ = awS«(*)«(-*)-6(*)6(-fc), if 9 > 0,6(|)^0 or g<0,a(-|) = 



and respectively 



^ = [ |, if 9 < 0,a(-f ) ^ or q > 0,6(f ) = 



3- if g > 0, 6(f) / or g < 0, a(-f ) = 



then the requirements above are satisfied for a e (k) while the symmetry condition takes the form 

1 n n ,a e (k)k-i/3 
r e (-k) = r e (k)==———. 3.16 

a e (k) k + i/3 

The symmetry (|3,16p yields a certain k — > —k symmetry of the solution M(x,t, k) of the RH 
problem (|3.6p . The relevant symmetry relation has been established in [5] in the case t = and 
x £ R . We shall need a version of that relation for the "complimentary" case, i.e. x = 0, and 
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t > 0. We shall perform the derivation in this case following practically the same arguments as in 
[5] (cf. the proof of Proposition 4.28 of [5]). 

Denote J(x,t, k) the jump matrix of problem (|3.6p . i.e., 



J{ X ,t,k) = e -«kx+2kh)* 3 j Q ( k y(kx+2k2t)c 



From (I3.16P it follows that 



J(0,t,-k) 
where 



WW 



-Oe(fc) fc+i/3„-4ifc 2 t 



1 



v r e(*0 al(k) fc-i/3' 

' (a e {k) \ [k-iP 



C-(k)J(Q,t,k)C+ l (k), (3.17) 



C(fc) = < 



a e (k) , 



+ 



a e (k) 



\ (k-ip 



A a e (k)J V * + 



o"i, Im k > 0, 



<7i, ImA; < 0, 



(3.18) 



1 

with o"i = ( ). This implies that the function 



(3.19) 



M(i, fc) := f (t, fc)M(0, t, -k)C{k) 

satisfies the same jump condition as M(0, Ac) does: 

M+{t,k) = M-(t,k)J(0,t,k), keR. 

The factor £(t,k) can be chosen so that M(t,k) — > J as — > oo and that M(t,k) has no 
singularities at = ±i/3. Indeed, introducing 







(3.20) 



with P(i) to be defined, one has 
Af(t,jfe) = (TiP 



fc+iff 



/ fp- x M(0, *,-&)) (P _1 M(0,t,-fc) 



' n 
A; -iff 

21 fc + i/3 



p- 1 M(0,t,-A:) 



12 I .D(/,)a,. 



22 



where 



(fa e (k) 



L>(fc) = diag{d 1 (A;),d 2 (A;)}= < 







i 



I V a e (k) / 
Therefore, there are no singularities at k = ±i/3 provided 



Im k > 0, 

— \ 

a e( fc ) I , ImA; < 0. 



(3.21) 



(p _1 M(0,t,i/3)) = and (p _1 M(0, t, = 0. 
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This suggests to determine P 1 (t) as follows: 



which gives 



with 



M 22 (0,t,i/3) -Mi 2 (0,t,i/3)\ 
M 2 i(0,t,-i/3) M n (0,t,-ip)) ' 

- 1 (M n (0, t, M 12 (0,t,i/3)\ 

1 j A(t) VM 2 i(0,t,-i/3) M 22 (0,t,i/3)y 1 ' ] 

A(t) = Mn (0, t, -i/3)M 22 (0, t, i0) - M12 (0, t , i/3)M 2 i (0, f , -i/3) 
= |Mn(0,t,-i/3)| 2 + |M 2 i(0,t,-i/3)| 2 >0. 

If the setting of the RH problem includes the poles, then, similar to the case considered in 
[S], one can verify that M(t,k) defined by (|3.19p - (|3.22p satisfies the same residue conditions as 
M(0, t, k) and hence the uniqueness of the solution of the RH problem gives 

M(t, k) = M(0,t,k), 

which reads in terms of M(0, t, k) only as 

M(0,t,-k) = a 1 P(t) ^]p-Ht)M(0,t,k)(^~ iP k l^D(k)a 1 . (3.23) 

Remark 3. The above arguments allow actually to prove the general symmetry formula which is 
valid for all x and t, 

M{-x,t-k) = a x P{t) ^]p- l (t)M(x,t,k)^~ iP k l^D{k)^ 

In the case t = and this is (up to the notations) the formula proven in [5]. 

We shall now show how the symmetry relation (|3.23[) alone can be used to establish the Robin 
boundary condition for u(x, t). To this end we first evaluate, using (|3.23[) . the entries Mu(0, t, —i/3) 
and M 2 i(0,t, -i/3). We have: 



Mn(0,t,- l/3 )= fc hm U) X) P ~ Hm ° ,t,k) ( k ^ k + v) m 



22 



P22(t) (p- 1 (t)M(0,t,-i/3)) 22 d 2 (-i/3) = ^M 22 (0,t,i/3)d 2 (-i/3) (3-24) 
-M-iP) 



Mn(0,t,-iy3)- 



A(t) ' 

where we have used the basic symmetry (|2.8|) . Similarly, we have 
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Pi 2 (t) (p- 1 (t)M(0,t,-i/3)) 22 d 2 (-i/3) = ^Mi 2 (0,t,i/3)d 2 (-i/3) ( 3 - 25 ) 
-da(-ij9) 



-M 21 (0,i,-i/3)- A(/ , 
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Comparing (|3.24p and (|3.25[) gives 



Mu(0,i, -i/3)(l-6>) = and M 2 i(0, i, — i/5) (1 + 0) = 0, (3.26) 

where 

= (3.27) 

which implies that either M n (0, i, -i/3) = or M 2 i(0, t, -i/3) = for all t > 0. 

The next (and the last) step is to explore an idea that was first suggested by A. Bobenko in the 
late 1980s, and which was used then in several works dealing with the algebro-geometric solutions 
of integrable equations (see e.g. [2]). 

Recall that V(t,k) := (Mn(0, t, k)e- 2ik ' 2t , M 2 i(0, t, k)e~ 2ik2t ) T satisfies the differential equation 
(|2.4j) with u = u(0,t) and u x = u x (0,t), i.e. 

+ 2iA; 2 ^i = i|u(0,i)| 2 #i + (2fcu(0,i) +iu x (0,t))¥ a> 

- 2ik 2 V 2 = -i|n(0,t)| 2 f 2 + (-2fcu(0,t) +iu x (0, 

From (13381 it follows that if *i(t, -i/3) = for all t > then -2i/3u(0,t) + iu x (0,i) = 0, or 
u s (0,i) -2/3u(0,t) = 0; if ¥ 2 (*,-i0) = then 2i,8u(0,t) +ro a .(0,t) = 0, or u x {0, t) + 2/3tt(0, t) = 0. 
Observe that, according to ([3. 15j) . /3 can be either g/2 or —q/2. But since the initial data satisfy 
the boundary condition ^(0,0) + qu(0,0) = 0, by continuity it follows that this condition holds 
for all t. 

A closer look at ()3.26j) and (|3.27p reveals that one can specify precisely whether (a) Mn(0, t, —i/3) 
or (b) M2i(0, t, —i/3) = occurs, depending on the sign of q and the properties of a(\\q\/2) and 
6(i|g[/2). Indeed, since A = \Mu{-ifi)\ 2 + [M 2 i(-i/3)| 2 > 0, the choice between (a) and (b) is 
determined by the sign of d 2 {— i/3). According to (I3.14j) and (13.15p . one can distinguish four cases. 

(i) If q > and b(iq/2) = 0, then /3 = § > and thus (see (EOT]! ) d 2 (-i/3) = a e (|). In turn, 

from (|3.14p it follows that in this case, a e (^) = o(^) > an d thus 1 + d2 (~ 1/3 ) > 0, which implies 
( see (HEMP ) that M 2 i(0, t, -i/3) = 0. 

(ii) If q < and o(-ig/2) ^ 0, then /3 = | < and thus d 2 (-i/3) = (a e (-| 
0. Hence, in this case one also has 1 + — A ^ > and thus M 2 i(0, t, —i/3) = 0. 



2 

> 



(iii) If 9 < 0anda(-ig/2) = 0, then /3 = -| > and thus d 2 (-i/3) = a e (-f ) = - 6(-f 



2 



< 0. 



-2 



Hence, in this case 1 — d2 > 0, which implies that Mu(0,t, —i/3) = 0. 

(iv) If q > and 6(ig/2) ^ 0, then /3 = -| < and thus d 2 (-i/3) = (a e (|)) 1 = - 6(f ) 

0. Hence, in this case one also has 1 - > 0) which implies that M n (0,t, -i/3) = 0. 

Summarizing, we see that M2i(0, t, — i/3) = corresponds to /3 = | while Mu(0, i, — i/3) = 
corresponds to /3 = — |, which is indeed consistent with the fact that (|3.28p implies u x (0,t) + 
qu(0, t) = for all t. 

4 Concluding remarks 

1. In the general Riemann-Hilbert approach to initial boundary value problems for integrable 
nonlinear equations [T2| [H [15] , an important step is the verification that the solution of the 
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underlying nonlinear equation obtained from the solution of the associated RH problem indeed 
satisfies the prescribed boundary conditions. In the general case, this can be done by mapping the 
master RH problem, in which the space parameter is taken to correspond to the boundary, to the 
RH problem associated with the t-equation of the Lax pair with a "potential" constructed from the 
prescribed boundary values, and showing that they are equivalent, in the sense that they produce 
the same "potentials". But this means that the latter RH problem must be well-defined, which in 
particular requires, in the case of semi-infinite time interval < t < oo, a precise description of 
the large-time behavior of the boundary values. On the other hand, such description is, generally, 
not available in full from the boundary conditions of a well-posed initial boundary value problem, 
which forces one to make certain a priori assumptions about this behavior. 

On the other hand, for linearizable boundary conditions, as we have shown on the example 
of the Robin boundary condition for the NLS equation, it is possible to verify directly that these 
boundary conditions hold, by using additional symmetry properties of an appropriately deformed 
master RH problem thus avoiding restricting a priori assumptions. The importance of this fact 
can be illustrated by the following simple observation: in the case q > 0, if b(iq/2) ^ (generic 
case!), then, as it follows from (I3.14p . a e (iq/2) = — |6(ig/2)| 2 < 0. On the other hand, a e (i£) — > 1 as 
£ — > +oo. Noticing that a e (i£) G R for all £ > 0, we conclude that a e (k) must have at least one zero 
for k G iM + . Consequently, as it follows from the asymptotic analysis similar to the one done in 
[5], stationary solitons generated by these zeros dominate the large-t asymptotics in the direction 
along the t-axis, which prevents from assuming the decaying behavior of u(0,t) and u x (0,t) needed 
for defining the associated spectral functions A(k) and B(k). 

2. For the initial data that do not produce zeros of a e (i£) and thus stationary solitons, the 
boundary value of the solution of the Robin problem decays as t" 1//2 , as it again follows form 
the asymptotic analysis similar to [5]. This is formally not enough in order to proceed with the 
approach of |10| . |15| . However, one can guarantee the needed rate of decay of the boundary value 
assuming that certain combination of the spectral functions a(k) and b(k), corresponding to the 
initial data uq{x) in the 1ST formalism, vanishes at k = to a high order. 

3. In our analysis of the IBV problem formulated for the domain x > 0, t > 0, we do not go 
"beyond the domain", working with x and t, as parameters in the RH problem, that stay in the 
domain prescribed by the problem. This is in contrast with the approach based on an appropriate 
continuation to x < with the help of the Backhand transformation allowing to control the necessary 
conditions at x = for all t [H \18\ [T9l [5] . Actually, folding the contour of the RH problem down 
to the real axis establishes the relationship between these approaches. Indeed, the functions r e (k) 
and a e (k) are exactly the reflection coefficient and the inverse transmission coefficient for the RHP 
problem on the whole axis associated with the initial value problem on the whole line with the 
Backlund-continued initial data, see |5j. 

4. An alternative way to see that u(x, t) constructed from the solution of the RH problem 
satisfies the boundary condition is based on using the further terms in the expansion of M(x, t, k) 
as k — > oo giving not only u(x,t) (as in (j2.16|) ) but also u x (x,t) [15J: 



u x (x, t) = lim A(k 2 M) 12 (x, t, k) + 2iu(x, t)(kM) 22 (x, t, k) 



(4.1) 



Indeed, let us substitute the expansion (for x = 0) 



M(Q,t,k) =1 + 



+ 



m 2 (t) 



+ ... 



k — > oo 



A; 
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into the symmetry relation (|3.23p taking into account that, in view of (|3.14|) and (j3. 15|) . 

a e (k) 



l + 0(£), if/3 = f, 



i-f + o(4), if/3 



Equating the terms of order k 1 one gets 











-2m\ 2 



On the other hand, from (I3.22D we have 

p a p-i _ 1 fMu(-ip) MuW 



A VM 2 i(-i/3) M 22 (i/3) 



0"3 



M 22 (i/3) -M 12 (i/3) 
-M 2 i(-i/3) M n (-i/3) 



(4.2) 



1_ (8 
A \0 -8 



where <5 = Mn(-i/3)M 22 (i/3) + M 2 i(-i/3)Mi 2 (i/3); here the fact that either Mn(-i/3) = or 
A^2i ( — = implies the diagonal structure of the resulting matrix. 

Further, in the case Mn(— i/3) = (recall that in this case (3 = — | ) one has 



5 = M 21 (-i/?)M 12 (i/3) = - M 21 (-i/3) 



M 2i (-i/3) 



and thus Pa^P 1 = — as- In the case M 2 i(— i/3) = (/3 = |) one has 



A 



Mu (-i/3) 



5 = M n (-i/3)M 22 (i/3) = Mn(-i/3) 

and thus Pa^P^ 1 = 03. In both cases, from fj4.2j) we have 

=m^(t) = 0. 

Now, equating the terms of order k~ 2 , we have for the off-diagonal part: 



(4.3) 





-2m 2 12 



-2m\\ 




which, in view of (|4.3p . reads 





m\ 2 {-2m\ x + iq) 



m\ 2 (t) = m\ 2 (t) ( -- 



in 21 (-2ml 2 - iq) 







(4.4) 



Finally, flUE) and (pHoT) . in view of and yield 

u x (Q,t) = Am\ 2 {t) = -2\qm\ 2 {t) 
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